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ABSTRACT 

We quantify the performance of a non-linear Wiener filter, constructed in wavelet 
space, at recovering some of the Fisher information that was lost in the weak lensing 
convergence field. The proposed method consists in a separation of the original field 
into the sum of a Gaussian and a non-Gaussian contribution. After filtering an en- 
semble of such fields, which are obtained from TV-body simulations, we find that we 
can recapture about four times more Fisher information, an effect that can potentially 
improve by a significant amount the constraining power of weak lensing surveys on 
cosmological parameters, including the dark energy equation of state uj. We compare 
this performance with that of a logarithmic mapping and find that the wavelet method 
can recover up to three times more information. 
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20ld : lHutereij|201(]| ). Because the light to mass bias is rather 
uncertain, the only direct procedure to weigh the universe 
is by measuring the deflection of light caused by the 
presence of matter between the source and the observer. 
In particular, the statistics of gravitational lensing can 
serve as a powerful probe of t h e mas s distribution of the 
universe (|Va n Waerta eke et all l200ll : IVafaei et all l2010l : 
iBeni amin et al . 2007). 

It was recently realized that weak lensing could also 
provide insight on dark energy via the measurement of the 
grow th function l|Huterei<200a : lAlbrecht et al.ll2006l : lHutereij 
l201[l[ ). and an international effort was put into motion in or- 
der to measure weak l ensing signal with unprecedent ac- 
curacy and resolution (ILSST Science Collaborations et al.l 



factor from clusters ( VoitI 2005ll and the detection of 
weak lensing signal (|Hoekstra et al.ll2006l ) have already set 
tight constraints on the dark energy equation of state lj 
(jOguri et al.ll200^ . and the weak lensing contribution pro- 
vides a complimentary approach as it is sensitive to dif- 
ferent systematic uncertainties. The goal of the combined 
dark energy experiments is to maximize a collective "flgure- 
of-merit", which influen ced he design of most of the future 
exper iments like L SSlt^ JLSST Science Collaborations et al. | 
20091) . EUCLIE0 ('Beau lieu et al.l |2010|V JDElVl^l (ICehrel 
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20061). SKAFI (Schihzz: 



I Schlegel et al.l 1200! 



20101) . CHIM EQ ( Pete rson et al.l 
120071 : iDewdncv et |2009| ). BOS: 
and PAN STARR^ 

It was soon realized in the BAO community that 
the constraining strength of these surveys depends di- 
rectly on the amoun t of Fisher information (|Fisheij Il935l : 
iTegmark et al.|[l997l ). i.e. statistically independent Fourier 
modes contained in the measurements of the power spec- 
trum bands, and one needs to maximize that information in 
order to minimize the uncertainty on cosmological parame- 



l2009l : fBeaulieu et al.ir2010l : lGehreIsll2010l ). Alternatives tech- 
niques su ch as the redshift dist ance measurements of su- 
pernovae ([Perlmutter et al] ll999'). the detecti on of baryonic 
acoustic oscillations (BAO) in galaxy surveys (jPercival et al.l 
l2007l : lE"isenstein et al.ll2005l ). the measurement of the growth 
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ters. Counting these Fourier modes is straight forward when 
the underlying density field is Gaussian, however departures 
from Gaussianity are to be expected, as arising f rom the 
non-l inear gravitational collapse of the density field (jHutereil 
|2002| ). In the theory of structure formation, large scales 
structures grow from an initially linear Gaussian random 
field, which progressively becomes non-linear through gravi- 
tational instabilities, starting from the smallest scales. Only 
the largest scales of the field remains intrinsically Gaussian, 
while the non-linea r Fourier modes start to couple together 
ijZhang et al.l |2003l). 

Rimes fc HamiltonI (|2005l . 120061 ) first measured the 
amount of Fisher information about the power spectrum am- 
plitude contained in the matter field, as a function of scale, 
from an ensemble of 400 A^'-body simulations. They found 
that in the largest scales, the Fisher information grows in a 
manner consistent with a Gaussian random field. However, 
they observed departures from Gaussianity in the trans- 
linear regime, in the form of an information plateau, followed 
by a second rise on much smaller scales. This was later in- 
terpreted in terms of the halo model as a transition between 
the two-halo and the one- halo terms ((Nevrinck et al.ll200d : 
iNevrinck fc Szapudill2007l ). 

This loss of Fisher information is an undesired effect, 
in the sense that it is equivalent to a reduction of the 
survey effective volume, and many strategies have been 
prop osed to re c over some of the erased cosmic informa- 
tion. IWeinberd (|l992l ) used a method called Gaussianiza- 
tion, which is a monotonic transformation of the smoothed 
galaxy distribution that reconstructs primordial density 
fiuctuations. Running A^-body simulations backwards in 



time l|Goldberg fc Sperge|200oi). or density field reconstruc- 
tion from linear theory ( Eisenstein et al.l 2007 _Noh et al] 



l2009l : iPadmanabhan et al.ll2009l : iNgan et al.ll201ll ') has also 
been successful at recovering parts of the lost information. 
More recently, it was found that a logarithmic transfor- 
mation of de nsity fields appears t o be effective on trans- 
linear scales (|Nevrinck et all |2009| ). INevrinck et"al] (|201ll ') 
have successfully recovered some of the Fisher information 
with a Gaussianization method that takes the Poisson noise 
into account. Non-linear Wiener filters, which can be de- 
signed to decompose a density field into Gaussian and non- 
Gaussi an parts, are also among the best techniques found 
so far ijZhang ct al. 2011). All of these techniques somehow 
attempt to diagonalize the covariance matrix of the matter 
power spectra by bringing back to the 2-point function some 
of the information that had leaked to higher order terms. 

The next step was to measure the impact of these non- 
Gaussianities on cosmological parameters. It was first shown 
that their inclusion has only a minor impact on the con- 
straining power about the BAO dilation scale, when both 
the power spectrum and the covariance matrix are obtained 
from A'-body simulations (Takahashi et al. 2011). In cur- 
rent data analyses, however, the power spectrum is often 
estimated wi th techniques that a s sume Gaussianity in th e 
matter field (|Feldman et all 1 19941 : IVogelev fc Szalavlll99d ). 
which have the unfortuna te effect of producing sub-optimal 
estimates of the mean (|Tegmark et al . 2006)- In the case 
where the survey select ion function is complex, the e stimate 
is likely to be biased (|Harnois- Deraps fc PenlbOllT ). In re- 
gards with the sub-optimal measurement, it was recently 
shown that the error bars on the BAO dilation scale that 



are consistent with a sub-optimal measurement of the power 
spectrum might be significantly larger, compar ed to those 
obtained under standard Gaussian prescription (|Ngan et al.l 
Soil). In the era of precision cosmology, these few percent 
level effects need to be considered when constraining dark 
energy parameters. 

In the pursuit of robustness and accuracy in weak lens- 
ing analyses, equal considerations must be granted to non- 
Gaussianities. Shear and convergence maps are indeed sensi- 
tive to the non-linear regime, typically at low and intermedi- 
ate distances. For instance, it was recently shown that Fisher 
information - about the amplitude of the lensing power spec- 
trum - contained in convergence fields w as also depart in; 
from the infor mation of a Gaussian field (jPore et al.ll200' 
ILu et akl lioiol). It was furthermore, shown that these devi 
ations could potentially impact the constraining power on 
the dark energy equation of state. The covariance matrix 
in the weak lensing angular power spectrum indeed shows 
strong correlation s across different scales, as confirmed by 
(|Seo et al.ll2011al ). 

Since then, much efforts have been made to Gaussian- 
ize the lensing fields as well, in an attempt to recover some 
information and thus improve the precision on c urrent and 
future measurements of cosmological parameters. (ISeo et al.l 
l2011al ) have also shown that the method of logarithmic 
mapping is able to pump back some of the information 
lost. Such transformation was shown to suppress the non- 
Gaussian contribution in the statistics of higher order cu- 
mula n ts, and can also supp ress the bispectrum ||Yu et al.l 
l201lf ). IJoachimi et al] l|201ll ) also showed that a Gox-Box 
transformation can also Gaussianize the fields and restore 
about the same amount of information lost in lensing fields, 
compared to logarithmic mapping. Those techniques basi- 
cally reconstruct a probability distribution function (PDF) 
of the K field which is much closer to that of a Gaussian 
jSeo et al.ll2011al : IJoachimi et al.ll201ll ). 

In the effort to Gaussianize the fields, non-linear Wiener 
filters are also a promising technique and work especially 
well with wavelet transforms whic h are well sui t ed fo r 
extracting multiscale information (|Fang fc Thewj Il998l ). 
Moreover, they were found to offer better performances over 
the s tandard Fo urier basis if the data are intermittent in na- 
ture (|Penlll999l ). They are proved to be successful in recov- 
ering Fisher information about the amplitude of the power 
spectrum of the dark matter field fehang et al.|[201ll ). and 
could potentially outperform other Gaussianization tech- 
niques of weak lensing maps. In this paper, we t hus construct 
a non -linear filtering method similar to that of IZhang et al.l 
(|201lf ) to Gaussianize the k fields, and we compare the in- 
crease the Fisher information contained in simulated maps 
with other methods. 

The outline of the paper is as follows: we first describe 
our strategy to construct weak lensing maps from TV-body 
simulations in We discuss in [jS] the application of the 
discrete wavelet transform (DWT) and Bayesian theory to 
filter out the non-Gaussian component of these fields. In i|4j 
we calculate the information content about the power spec- 
tra in the Gaussianized k fields. Discussion and conclusion 
are presented in ^ 
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Figure 1. Random selected k field constructed from A^-body simulation (left panel). After the non-linear wavelet Wiener filtering, this 
field is decomposed into Kg (middle panel) and Kng (right panel). 



2 SIMULATIONS AND CONSTRUCTION OF 
CONVERGENCE MAPS 

We ran a series of 185 A''-body simulations to ge nerate con- 
verge nce maps, using the fast parallel cubepSm (|Merz et al.l 
[2003), a particle-mesh Poisson solver that has sub- grid res- 
olution, thanks to the calculation of the exact Newton force 
between particles. It was run with 1024^ particles on eight 
IBM nodes of the Tightly Coupled System on SciNet, each 
node being equipped with 32 cores, 12 8Gb of RAM, wit h 
an infiniband connection across nodes l|Loken et al.ll201oh . 
We output a series of periodic surface density on 2048^ pixel 
maps, along the three orthogonal directions of the cube at 
each specified redshifts. 

The simulations were designed to optimize the usage of 
the simulated comoving volume in the construction of the 
past light cone, which consist of a juxtaposition of cubes of 
231.1Mpc//i per side for z > 1, and of 147.0Mpc//i per side 
for lower redshifts. The hght cone opening angle was set to 
3.58 degrees, which exactly touches the edges of the small 
box 'At z = 1, then propagates into the larger boxes until 
z = 2, beyond which we used the periodicity of the simula- 
tions to populate the volumes. This creates repeated struc- 
tures, which increase the correlation across different Fourier 
modes in the mass density at the percent level. But this 
has negligible impact on our results since the lensing kernel 
strongly suppresses the contributions from such high red- 
shifts. At most, such increased correlation indeed accentu- 
ates the non-Gaussian features, but this effect gets strongly 
suppressed by the projection along the line of sight. 

Simulations started at an initial redshift Zi — 40 for the 
lower redshift boxes and 200 otherwise. The cosmological 
parameters used are flu = 0.279, Qa = 0.721, with Hub- 
ble constant h = 0.701, and we obtained a transfer function 
from CAMB with Us = 0.96. The power spectrum normaliza- 
tion was then specified to be as = 0.815. 

The convergence k field is obtained from the projection 
of the matter overdensity S along the line of sight 9, weighted 
by the lensing geometry and, potentially, a source-galaxy 



distribution. It can be expressed as 

^io,x)= rw{x')5{x',r{x')e)dx', (1) 

Jo 

where x is the comoving distance in unit of c/Ho, and Ho — 
100 h km s"^ Mpc"^ The weight function W{x) is 

W{x) = ^^lMg{x){l + ^) (2) 

where g{x) is determined by the source galaxy distribution 
function n{z) and the lensing geometry: 

9ix) = r{x) n{x!)'^^^Ax! ■ (3) 

Here, r(x) is the radial coordinate and is equal to x for the 
fiat geometry we consider. 

The density fields are converted into convergence maps 
by stacking the images with the appropriate geometrical 
weights W{z), through the comoving volume contained in 
the pa st hgh t cone. This tiling method was developed in 
iSeliakI l| 19981 ) and assumes both the thin lens and Born ap- 
proximations. At each lens plane, we choose randomly one 
of the y- and z- directions, and we shift the centre of the 
plane. This effectively suppresses much of the correlation 
that exists across the lenses. We then interpolate the lens 
onto a pixel map of constant angular size. Since we chose a 
uniform galaxy distribution, and a source placed at Zs — 3.0, 
most of the lensing contribution comes from 2 ~ 1.5. One of 
these K map is shown in the left panel of Fig[T] 



3 WAVELET NON-LINEAR WIENER 
FILTERING OF k FIELDS 

In this section, we first briefiy review the properties of the 
wavelet transform that are relevant to our discussion, then 
we describe our non-linear Wiener filter, and finally we 
quickly present the log transform, which has been previ- 
ously used to Gaussianize the k fields and against which we 
compare our results in section |4] 
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3.1 Discrete Wavelet Transform (DWT) 

Similar to the Fast Fourier Transform (FFT), the Discrete 
Wavelet Transform (DWT) is an invertible, linear opera- 
tion, and can be considered as a rotation in function space. 
We consider in t his paper t he Daubechies-4 (hereafter DB- 
4) wavelet basis |Daubechies 1992), which contains families 
of scaling functions </> and difference functions (or wavelet 
functions) ip that are orthogonal, continuous and have com- 
pact support. In a DWT decomposition, the data which, in 
our case, have 2^ grid elements per dimension, are thus ex- 
panded into a combinations of these orthogonal basis, and 
weighted by wavelet coefficients e. 

In two dimensions, each wavelet coefficient thus depend 
on two scale indices [31,32) - controlling the dilation of the 
wavelet DB-4 functions - and two location indices (lijh) - 
controlling its translation. On a given dimension, the grid 
scale corresponding to a specified dilation is L/2-' , where 
L = 2048 in our case. Similarly, the index / corresponds to 
the spatial location comprised in the range IL/2-' < x < 
{I + l)/2-' . The 2-dimensional convergence field k{6{xi, X2)) 
can be expanded as 

1 1 

k(xi,X2) = ^ ^ eo,0;ii,i2</'o,ii (a;i)0o,i2(a;2) 



+ J2^fl,0;h,l2'l'jlMixi)<l'0,h{x2) 
jl=0 ii=0 i2=0 
1 J-1 2^2 -1 

+ Z] 5Z 5Z 4j2;ii.i2''^o.'i(^i)^.'2,i2(a^2) 

Il=0j2=0 l2=0 
J-1 2J1-1 J-1 232-1 

- 5Z 5Z X] 5Z ^Ji.J2;ii,i2V'ji,ii(2;i)fe,i2(a;2), 

J1=0 li=0 J2=0 1-2^0 



where 



■PjA^) = \Ij^^{2'x/L~1) 



t^,Ax)^\j-^H2'x/L-l). 



(4) 



(5) 



(6) 



With the combination of these two basis functions, scaling 
function coefficients (hereafter SFC's) e 's and three kinds 
of wavelet function coefficients (hereafter WFC's) e 's can 
be calculated by 

<^ji.32;h.h ^ JJ f^{xi,X2)(t)ji,h{xi)(j}j2,i2{x2)dxidx2, (7) 
^fi,j2;h,i2 ^ ff f^{xi,X2)i'ji,h{xi)(p32.h{x2)dxidx2, (8) 
fi,j2;h,i2 ^ ff fiixi,X2)(l)j,,i,{xi)^j^,iAx2)dxidx2, (9) 



and 



^ii,j2;h,i2 = J J Plixi,X2)i^n.ii{xi)'ipj^,i^{x2)dxidx2. (10) 

For each simulation, the n field is thus wavelet trans- 
formed, and each of the four kinds of coefficients found in 
Eas.((7l fT0|) are stored in a 2-dimensional field, preserving 



the gr id resolution fsee iFang fc Thew3 l| 19981 '): [Press et all 
(|l992l ) for more details). 



3.2 Non-linear Wiener Filtering 

Our strategy to construct a non-linear Wiener filter relies 
on the fact that in wavelet basis, the non-Gaussianities are 
clearly characterized in the PDF of the WFCs Cj^j^'.h.hy 
which we obtained from EQ. pO|) . We thus construct our filter 
by splitting the wavelet transform of the original map, which 
we label K, into two components: a Gaussian (G) and a non- 
Gaussian (TV) map. Namely, in wavelet space, we have 



K = G + N. 



(11) 



Since wavelet transforms are linear operations, we can in- 
verse wavelet transform the above equation and write, in 
real space. 



(12) 



where the original map (k) is expressed as the sum 
over a Gaussian contribution (hereafter Kg) and a non- 
Gaussianized contribution (hereafter Kng). Our goal is thus 
to design a filter that concentrates most of the collapsed 
structure in Kng, and thus produces Kg that are closer to 
linear theory. We perform this operation on our simulated 
maps, compute their power spectrum, construct a covariance 
matrix and measure the Fisher information of both compo- 
nents separately. Then we can finally compare our results 
with the unfiltered maps. 

The algorithm we use to devise the filter is a 2 - 
dimensional version of that presented in (|Zhang et al.ll201ll ) , 
which we briefiy describe here again for completeness. The 
filter acts on each wavelet mode, which is defined as a collec- 
tion of all WFCs having the same two scale indices (ii,j2). 
Let us focus on a given wavelet mode, or a pair of scale 
components, or dilations (ji,j2), whose value at location 
(^1,^2) is labeled x for simplicity. We construct the PDF of 
the WFCs in this wavelet mode, f{k), by looping over the 
other two indices {h, h), and express this measurement as a 
convolution of the PDF of two components. The Gaussian 
contribution G(x) = K{x) — N{x) has a Gaussian PDF by 
assumption. We label the PDF of the non-Gaussian contri- 
bution as 0{n), and write 



f{k) 



0(n)exp( {k — n) )dn 



(13) 



This description also assumes that all the measurements of 
X are statistically independent. Thanks to Hayes's theorem, 
we can then calculate the conditional probability P[N\K) — 
P{K\N)P{N)/P{K). The posterior conditional expectation 
value is written as 

{N\K = fc) = I exp[-i(n - k)^]e{n)ndn 



= K+ ^= dk 

V2^f{k) 

= K+ {In f )'{k) 



exp[— — (n — fc) ]0(n)dn 



(14) 



The expectation value of the Gaussian distribution is simply 
recovered from Ea. (|lip . In other words, knowing f{k) allows 
one to solve for both G and A'^ given K. 

When measuring the PDF of a given wavelet mode 
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Figure 2. WFC PDFs (red dots) and their analytical fits by Eq. I|15|I (black solid lines). Their corresponding Gaussian PDFs, which 
have the same standard deviations with the central regions of the fitted PDFs, are also shown in blue dashed lines. The upper left, 
upper right, bottom left and bottom right panels show four different wavelet modes, from smaller (more non-Gaussian) to larger scales 
(more Gaussian), as indicated by the value of (ii,i2)'s in the titles. They are roughly corresponding to multipole moment £ ~ 9.1 X 10^, 
2.3 X 10"^, 5.7 X 10^ and 2.8 X 10^, respectively. The two parameters a and s for Eg.l llSI I are also shown in the title of each panel. 



(jitji) in practice, we encounter the difficulty of poorly de- 
termined PDFs, and thus even worse are the filter functions. 
This condition is rooted in the scarceness of the WFCs (espe- 
cially for larger scale wavelet modes, containing less WFCs) 
and in the instability of numerical differentiation, fn order to 
resolve this issue, we construct f{k) for each wavelet mode 
by looping over the spacial indices [h, I2), and also over all 
the realizations. We further improve the resolution of the 
PDF of each mode with a two parameter analytical func- 
tion: 



/pDp(a::) 



r(ia.^ 



r(ias2-i) 



/ 2 



(15) 



which follow the PDF curves (see Fig[2]), where the F's are 
the usual Gamma functions. The two parameters a and s 



are actually extracted via their relationship to the second 
moment and the fourth moment 1714 of the PDF: 



5m4 



9mi 



2m2m4 



and 



2m2m4 



m4 



(16) 



(17) 



Having measured the expectation value of both compo- 
nents, we can now loop back over all spatial indices {h,l2) 
and separate each coefficient into the two components, by 
applying the following (non-linear) Wiener filters: 



«;g(x) = G/K = 



(ln,f)'(fc) 



(18) 
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^ng(x)^iV/j^=l + (l'^^J(^) (19) 

which are functions of K only, and, we recall, apply solely 
to the current wavelet mode (ji,j2)- 

With the above parametrization in terms of (a,s), the 
Gaussian filter function can be expressed by combining 
Eq.lTHI) with Eq.jrn}: 

w^{x)^{l + ^^ . (20) 

Note that the final filter function depends only on s, which 
characterizes the extent of the departure from a Gaussian 
PDF, i.e., the greater the s, the smaller departure from 
Gaussian. 

In the limit of a pure Gaussian PDF 

/^DF=AA(0,a2) = -^e-^^, (21) 

the second moment and the fourth moment are given by 
rni = o-^ and = 3o-*. From Ea. p6ll7|l we have a = cr"^ 
and s — >■ +00, which is understandable. In this case, from 
Ea. l|18p the filter function w'^ reduces to a constant unity 
function w'^{x) = 1. This is to be expected since for those 
Gaussian or nearly Gaussian distributed PDFs, we do not 
need to filter at all. 

In Fig[2] we select four wavelet modes, characterizing 
four different scales of the k fields, and plot the PDF and 
fitted PDF for each of their WFCs. In order to visualize 
better the non-Gaussianities of the PDF, we also plot for 
each panel a Gaussian PDF which fits the central region 
of the fitted PDF Eg. lfTSl) . As indicated by the dilations 
Oil jz) in each of the panel's title, the corresponding scale 
increases when one looks from the upper left panel to the 
bottom right panel. The corresponding multipole moments 
for those four chosen {ji,j2)'s are roughly f ~ 9.1 x 10"^, 
2.3 X 10^, 5.7 X 10^ and 2.8 x 10^, respectively. One can 
also see that, as the scale goes up (or £ goes down), the 
parameter s grows larger, and the PDF (black solid line) of 
that wavelet mode becomes more Gaussian, and closer to 
the Gaussian distribution (blue dashed line). 

Note that the wavelet non-linear Wiener filter is a 
parameter-free method. In each wavelet mode, the filter is 
determined only by the 1-point PDF of all the WFCs. The 
method is stable in that, if we make little changes in data in 
real space, then all the WFC's PDFs will have little change, 
but will not have much effect on the final fit, thus will not 
change the non-Gaussian decomposition. In contrast, if one 
adds in a lot of non- Gaussian features on certain scales, then 
the PDF of the WFCs, especially on that scale, will be more 
non-Gaussian. This results in a stronger filtering and those 
structures will be filtered out. 

We finally repeat this process for all wavelet modes, 
constructing the filters Wg and Wng and separate each « map 
into a Gaussian and non-Gaussian part. As an example, one 
decomposition among those 185 realizations is shown in the 
middle and right panel of Fig[T] 

3.3 Logarithmic Mapping of k Fields 

In order to quantify the performance of our method, we 
wish to compare the recovery of Fisher information in our 




Figure 3. Probability distribution function of k, Kg, Kng and Kin 
field. For comparison, A standard Gaussian distribution's PDF is 
also shown with a dotted line. 

Gaussianized fields with t hat of a compet ing method. We 
follow the prescription of (|Seo et al.ll2011al ). and "log-map" 
each of our k fields 

Km = Ko In (l + ^) . (22) 

where kq is defined so as to keep the argument of logarithm 
positive, i.e. kq < min(K), with min(«;) being the minimum 
value of the k field. This transformation is designed such that 
the non-linear peaks, who show strong convergence values, 
are attenuated. As a result, higher-order statistics become 
less important, which translates in a decrease of the non- 
Gaussian contribution. We can characterize the degree of 
alteration of the field by defining r = ko/| iiiin(K)|, where 
1 < r < 00. The smaller the r, the more the log-mapping 
alters the field. In order to draw general conclusions, we 
sample the r space and try different values in the following 
sections. These log-transforms are also applied onto wavelet 
filtered k fields: 

Kln+g = KO In ^1 + , (23) 

The results of both methods and of their combination are 
discussed in the next section. 



4 ANALYSIS 

In Fig[Tl we can see by eye that the peaks in the unfiltered k 
maps are mostly transferred into Kng, leaving Kg with much 
less structure. We also plot the real-space 1-point probabil- 
ity distribution function in Fig[31 The result shows that the 
original convergence map has a very large skewness. We ob- 
serve that the logarithmic mapping method have the best 
effect in recovering a Gaussian PDF, while the non-linear 
Wiener filtering mildly removes some of the skewness. This 
shows that the wavelet techniques is optimized for restoring 
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Figure 5. Cross-correlation coefficient matrices as found from 185 angular power spectra of k (left panel), Kg (middle panel), and Kng 
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4.1 Angular Power Spectra 

We calculate for each map their dimensionless angular power 
spectra C(t), which is defined as 



i{i+i)p{i) 

27r 



(24) 



where I is the multipole characterizes the scale in the two- 
dimensional K-fields. We construct 20 ^-bins and for each 
wavenumber \l\ , the power spectrum P{£) is estimated by av- 
eraging all the Fourier modes within the bin, while £ = (^bin 
is determined by averaging all \£\ that fall into the bin. In 
Fig[3]we plot the mean angular power spectra and error bars 
of K, Kg, and Kng fields respectively. We can see that on lin- 
ear scales, where £ < 10^, the Gaussianized power spectrum 
is nearly unchanged. On non-linear scales, however, it drops 
by a factor of two, while the error bars are reduced. Hence, in 
the presence of moderate noise, such a lowering of the power 
spectrum shall not affect significantly the detectability. 



4.2 Cumulative Information 

We measure the covariance matrix of the C{£), which cap- 
tures the correlation between the variance, or the error bars, 
of the power spectra at different scales £. If the measure- 
ments were completely uncorrelated, the diagonal of the co- 
variance matrix would be the variance at each value of £, and 
all the off-diagonal entries would be zero. Mathematically, 
the covariance matrix C is defined as 

1 ^ 

C{£,£') = ^^^[C,(£) - {Cmm') - (C(/))], (25) 



Gaussian features in the 2-point function, whereas logarith- 
mic mapping is more effective on the PDF. Both of these 
statistical estimators are useful in data analysis, as they are 
sensitive to different systematics and probing the cosmology 
in slightly different ways. 



where is the number of realizations and {C{£)) is the mean 
angular power spectrum over all realizations. Note that the 
covariance matrix C {£,£') is not to be confused with angular 
power spectrum C{£). 

The cross-correlation coefficient matrix, or for short the 
correlation matrix, is a normalized version of the covariance 
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matrix, where each value is divided by the square root of 
the diagonal values as follows: 



C{£,£') 



(26) 



The three correlation matrices for our unfiltered, Gaussian- 
ized and non-Gaussianized k fields are shown in the three 
panels in Fig[S] As expected, the original correlation matrix 
shows cross-correlation greater than 60% and up to 80% 
where £ > 10^, while on larger scales, which is still in the 
linear regime, the matrix is diagonal. In the Gaussianized 
correlation matrix, the correlated region is much smaller, 
so the correlations between different Fourier modes are sup- 
pressed to some extent. In contrast, in the non-Gaussianized 
matrix, the cross-correlation is even higher and even spread 
into linear region. 

The cumulative, or Fisher, information, is a measure- 
ment of the number of independent Fourier modes presented 
in a field up to a given scale, here £„. It is obtained as fol- 
lows: for a given wavenumber £„, we select the subsection 
of the covariance matrix up to that scale, then we invert 
this sub-matrix and sum over all its elements. The results 
are meaningful when using a normalized covariance matrix, 
which is defined as 



C'norm(^,^ ) — 



C{£,£') 



m{C{£'))' 
and the function of cumulative information is 



£. 



(27) 



(28) 



We plot the cumulative information contained in the un- 
filtered, Gaussianized and non-Gaussianized angular power 
spectra respectively in Fig[6l We also over-plot the informa- 
tion curve obtained from Kin with r — 1.5 which seems to 
have the best result among log transforms, and the theo- 
retical Gaussian information. In analogy with the Fisher in- 
forma t ion m e asured in the densi ty field iRimes fc HamiltonI 
l|2005l . l2006l ): IZhang et al] l|201ll ). K-fields' information also 
deviates from the Gaussian prediction towards £ = 800. 
However, in the weak lensing case, the departure from Gaus- 
sian predictions is somehow attenuated, an effect of the an- 
gular projection across multiple scales. 

Note that because we are using a finite number of sim- 
ulations to measure several bins of data, even for a Gaus- 
sian ra ndom field, the inve rse of the covariance matrix is 
biased l|Hartlap et al.l [20071 ). The biases are bin-dependent 
and are up to 10%, and would be applied to all the infor- 
mation curves, but we are interested only in the difference 
between those curves. For this reason we do not include this 
bias in all calculations. 

We also note that our unfiltered Fi sher in formation is 
somewhat lower than that presented by iDore'e t al. (2009), 
but this apparent discrepanc y is caused b y a di fferent choice 
of galaxy window function. iDore et al.l (|2009l ) opted for a 
series of top hat windows, extending from z ~ 1.0 to 2: ~ 3.0, 
while we include galaxy counts all the way down to z = 0. 
This effectively enhances the amount of non-linearities in 
our fields, thus reducing the Fisher information. 

After performing the wavelet non-linear Wiener filter- 
ing and logarithmic mapping method, both the Kg and Kin 
fields recovered information, however, the wavelet technique 



restores nearly 4 times more information than the logarith- 
mic mapping for £ > 10000. This performance gap gets even 
larger for r — 3.0, r — 2.0 or r = 1.05. We should also 
mention that wavelet method has the advantage to be a 
parameter-free method, in the sense that degrading the res- 
olution by factors of two does not affect the WFCs PDF es- 
timation of other scales. For the logarithmic method, how- 
ever, one needs to re-optimize the value of r for different 
resolution, which affects the Fisher information at all scales. 
For example, there are likely coarser resolutions where the 
log transform would perform better (and worse) than the 
resolution shown here. 

We also investigated whether these two methods could 
be combined together such that the total information recover 
exceeds that of the two methods separately and found that 
these two methods should not be used in conjunction. When 
we perform a logarithmic mapping after wavelet filtering, the 
recovered Fisher information is lower than for the wavelets 
alone. This is likely caused by the fact that a logarithmic 
mapping applied on a nearly Gaussian field transforms the 
field into something else, even less Gaussian. Because the 
wavelet transform is resource consuming, we have not tried 
the other way around, i.e. performing the logarithmic map- 
ping first, which would involve a filtering for each sampling 
of kq. However, we have strong reasons to believe that the 
gain would also be minimal, if not worse than the wavelet on 
its own. As discussed before, we observe that the PDFs of 
the log-transformed fields, in both real and wavelet spaces, 
are nearly Gaussian. This means that the wavelet filtering 
would have a very small impact, as most of the contribution 
would directly fall in the Gaussian filter. As a result, the in- 
formation that is recovered by the wavelet non-linear Wiener 
filter, but not by the log-transform, will not reappear with 
an additional Wiener filtering. However, other Gaussianiza- 
tion methods, especially those that attempt to go back in 
time - and thus only leave some non-Gaussian features - are 
likely to combine much better with the wavelet non-linear 
Wiener filters. For example, we are in the process of testing 
whether wavelet filtering can be used in conjunction with 
a density reconstruction algorithm to optimize the Fisher 
information. 

In addition, it has been shown that the logarithmic 
mapping has very little gain in Fisher info rmation when on e 
includes realistic levels of Gaussian noise (|Seo et al.ll2011bl ). 
which are inherent to weak lensing observations. Although 
we have not yet tested the performance of wavelet filtering 
in such a noisy environment, it has for sure the advantage 
of being able to extract a non-Gaussian component, which 
still contains most of the collapsed structures, and give us a 
novel handle in the data analysis. 

The extension of this work, which we postpone to a 
future paper, is to propagate the uncertainty of the mea- 
sured angular power spectrum onto that of cosmological 
parameters, fo l lowing the Fisher formalism summarized in 
lAlbrecht et al.l (|2006 ). and to compare the constraining per- 
formances of the covariance matrices from the unfiltered and 
from the Gaussianized fields to the simple Gaussian anal- 
ysis. Quantifying the difference in uncertainty on the dark 
energy figure-of-merit is the final objective, and a significant 
improvement would lead to a enhanced yet robust precision 
on uj. 
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5 CONCLUSION 

We ran a series of A'^-body simulations to generate weak lens- 
ing convergence maps. After analyzing their angular power 
spectra, as expected, we find that their cross correlation ma- 
trix is highly correlated on small scales {£ > 10^), since these 
scales are in non-linear regime. This non-linear clustering 
makes the Fisher information contained in the angular power 
spectra less than Gaussian information. The lost informa- 
tion can be obtained by higher statistics like bi-spectrum 
(three point correlation function), tri-spectrum etc., but to 
avoid dealing with tedious higher order calculations, we use 
two mathematical tools to Gaussianize k fields, namely a 
logarithmic mapping and non-linear wavelet Wiener filters. 

The method of wavelet non-linear Wiener filtering and 
logarithmic mapping can both increase the Fisher informa- 
tion in the angular power spectra, but wavelet method has 
better results with a gain of a factor 4 ai £ — 10000. How- 
ever, one needs to keep in mind that with the forecasted level 
of sh ot noise of fu ture survey lik e EUCLID (|Beaulieu et al.l 
I2OI0I ) and JDEM (|Gehrelsll20l"ol ). the relative importance of 
non-G aussianities will be reduced significantly (|Dore et al.l 
I2OO9I ). We also find that these two methods are not comple- 
mentary, and that their combination reduces the information 
recovery compared to a sole Gaussianization of the fields. 
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